Long-range Coulomb interactions in surface systems: a first principles description 
within self-consistently combined GW and dynamical mean field theory 
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Systems of adatoms on semiconductor surfaces display competing ground states and exotic spec- 
tral properties typical of two-dimensional correlated electron materials which are dominated by a 
complex interplay of spin and charge degrees of freedom. We report a fully ab initio derivation of 
low energy Hamiltonians for the adatom systems Si(lll):X, with X=Sn, Si, C, Pb, that we solve 
within self-consistent combined GW and dynamical mean field theory ("GW+DMFT"). Calculated 
photoemission spectra are in agreement with available experimental data. We rationalize experimen- 
tally observed tendencies from Mott physics towards charge-ordering along the series as resulting 
from substantial long-range interactions. 

PACS numbers: 71. 15. Mb, 73.20. At, 71.10.Fd, 71.30.h 
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Understanding the electronic properties of materials 
with strong electronic Coulomb correlations remains one 
of the biggest challenges of modern condensed matter 
physics. The interplay of derealization and interactions 
is not only at the origin of exotic ground states, but also 
determines the excitation spectra of correlated materi- 
als. The "standard model" of correlated fermions, the 
Hubbard model, in principle captures these phenomena. 
Yet, relating the model to the material on a microscopic 
footing remains a challenge. Even more importantly, the 
approximation of purely local Coulomb interactions can 
become severe in realistic materials, where long-range in- 
teractions and charge fluctuation physics cannot be ne- 
glected. 

Systems of adatoms on semiconducting surfaces, such 
as Si(lll):X with X=Sn, C, Si, Pb, have been suggested 
[H to be good candidates for observing low-dimensional 
correlated physics. Commonly considered to be realiza- 
tions of the one-band Hubbard model and toy systems 
for investigating many-body physics on the triangular 
lattice, such surfaces have been explored experimentally 
and theoretically fToM3ll ] .These so-called cv-phases 
show a remarkable variety of interesting physics includ- 
ing commensurate charge density wave (CDW) states 
@; H Hand isostructural metal to insulator transitions 
(MIT)[14|. However, while specific systems and/or phe- 
nomena have been investigated also theoretically, a com- 
prehensive understanding including materials trends is 
still lacking. A central goal of our work is to present a 
unified picture that relates, within a single framework, 
different materials (adatom systems), placing them in a 
common phase diagram. 

We derive low-energy effective Hamiltonians ab initio 
from a combined density functional and constrained ran- 
dom phase approximation (cRPA) scheme [32| in the im- 
plementation of [33| (see also the extension to surface sys- 



tems in [34[). While the first surprise are the relatively 
large values of the onsitc interactions which we find to be 
of the order of the bandwidth (« 1 eV), most importantly 
we show that non-local interactions arc large (nearest- 
neighbor interaction of ss 0.5 eV) and, hence, an essential 
part of the resulting many-body Hamiltonians. This re- 
sult confirms previous speculations about the importance 
of non-local effects in these materials [2ll [2^ . We solve 
these Hamiltonians within fully self-consistent combined 
GW and dynamical mean field theory ("GW+DMFT") 
[35l | . calculating in particular (single particle-) angular 
resolved photoemission spectra (ARPES) and the (two 
particle-) charge susceptibility. We identify a clear- 
cut materials trend starting from Si(lll):C deep in a 
Mott phase to Si(lll):Pb which shows tendencies to- 
wards metallicity and charge-ordered states driven by 
non-local interaction terms. Comparing our results to 
available experimental data yields encouraging insights: 
Without adjustable parameters we reproduce the experi- 
mentally measured gap size of insulating Si(lll):Sn and 
its transition to a "bad-insulator" at elevated tempera- 
tures. Moreover, based on the charge susceptibility, we 
identify the electronic tendency of Si(lll):Pb towards 
charge-ordering of the so-called 3x3 symmetry, which 
is indeed seen experimentally by scanning tunneling mi- 
croscopy. Our work is the first one that addresses the 
electronic properties of real materials on the basis of fully 
self-consistent GW+DMFT calculations (for a non-self- 
consistent calculation se e 13(1 , for self-consistent calcula- 
tions for models see HB-Hg). [13] 



The single-particle part of the Hamiltonian is calcu- 
lated in the local density approximation of density func- 
tional theory. In Fig. [1] we present LDA bandstructures 
for the series Si(lll):{C, Si, Sn, Pb}. For all considered 
systems the surface-state in the semiconducting gap is in- 
deed responsible for a well- separated, single band around 
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FIG. 1. (Color online) Bandstructures of the a- \/H x \/3 
phases for Si(lll):X with X=Sn, Si, C, Pb [H]. The color 
of the bands denotes their respective orbital character. Red 
color indicates a p z -like "apical" character, while the blue 
color denotes p x , H -like (i.e. planar) character. The black dots 
represent the tight binding fit given by Eq. |T]) and hopping 
parameters from Tab. U 



the Fermi energy. In red (gray) we plot the contribu- 
tions stemming from the p z -orbital of the adatom while 
we plot the adatom p X J ,-character in blue (dark gray). 
Even though the actual molecular orbital composition 
might be complicated, the half-filled surface band has a 
clear-cut "apical" (i.e. carrot-like) character. For our 
calculations presented below we directly use the ab ini- 
tio derived dispersion relation. However, for the purpose 
of analysis we note that the tight-binding dispersion of 
the half-filled surface band can be well fitted using up to 
third-nearest-neighbor hopping (t, t' , and t") by: 

£k = 2t ■ (cos(k x ) + 2 cos(fc. T /2) cos(v / 3/2fc y )^ 
+ 2t' ■ (cos(V3k y ) +2cos(3fc x /2) cos(V3/2fc y )) 

+ 2t" ■ (cos(2/c x ) + 2 cos(fc x ) cos(V3kyfj (1) 

The values for the hopping integrals can be found in 
Tab. U and we plot the analytically calculated bands in 
Fig. [1] as the black dashed line. The quality of the fit 
supports the picture of Wannier-like orbitals with a fast 
decaying real space overlap on neighboring sites. 

In order to determine the interaction parameters as 
partially screened matrix elements of the Coulomb in- 
teraction within the cRPA one has to choose a suitable 
energy-window around the Fermi energy encompassing 
the surface band. The bare interaction parameters are 
calculated by means of explicit evaluation of the ra- 
dial (Slater-) integrals of the Wannier functions. Sub- 
sequently, the dielectric tensor is obtained within cRPA 
for local and non-local interaction parameters |34| . The 
results are summarized in Tab|I] 

The bare onsite interaction parameters (Vq) vary be- 
tween 6.0 eV for Si(lll):C and 4.3 eV for Si(lll):Pb 
decreasing monotonously within the series. The onsite 
Uq is reduced roughly by a factor of 4 — 5 due to cRPA 
screening. At first glance the onsite/static Uq of the order 



TABLE I. Values of the bare (V) and static, screened (Uq = 
U(iu = 0)) values for on- and intersite nearest neighbor (nn) 
interaction parameters. Also reported are the values of the 
static component of the effective U(ui — 0) calculated from 
GW+DMFT, see text. 
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of 1 eV - about twice the size of the bandwidth - strongly 
points towards Mott physics. This is, however, a prema- 
ture conclusion due to the effect of non-local interaction 
terms. The first non-local contribution (nearest-neighbor 
interaction) U\ [bare V\] is 0.5 eV [2.8 eV]. Remarkably, 
the value is - opposed to Uq [Vb] - almost the same for all 
materials. The reason is that the intersite overlap of the 
orbitals is so small that the Coulomb energy corresponds 
to the electrostatic energy of two point charges. With the 
virial theorem (E tot -) = \/2(V), we quantify this argu- 
ment by a rescaled hydrogen problem with effective Bohr 
radius of 6A (^distance of adatom sites): 
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which roughly matches the value of our bare intersite 
interaction parameters. The second, likewise remark- 
able, observation is that the screened values U\ are ex- 
tremely close to the value we get by assuming a static 
continuum approximation on the surface of a dielectric 
medium: Vi/e^ 1 , where e|" rf = ^(esi + 1) is the static 
dielectric constant of silicon on the surface. The reason 
is straightforward: The adatom distance (6A) is already 
large enough compared to the atomic structure of the sil- 
icon substrate (« 2A) so that local field effects (included 
in cRPA) are negligible. Following this reasoning we can 
calculate longer range interaction terms by simply scal- 
ing U\ with a/r, i.e., with the distance in units of the 
nearest-neighbor distance a, i.e., U-2 = Ui/V3 and so on. 
In this respect Ui is not only the nearest-neighbor inter- 
action, but the parameter that quantifies the strength of 
non-local interaction. 

To solve the effective low-energy Hamiltonians result- 
ing from our parameter-free downfolding procedure we 
implement the combined GW+DMFT scheme [35|, 42 1 



and calculate spectral properties and charge-charge re- 
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FIG. 2. (Color online) Momentum-resolved spectral function 
at T = 116K of Si(lll):X with X=Sn, Si, C, Pb obtained by 
analytical continuation of GW+DMFT imaginary-time data. 
The Fermi energy is set to ep = and indicated by the white 
dashed line. On the bottom right we show the spectral weight 
at the Fermi energy as a function of temperature. 



sponse functions. Fully self-consistent GW+DMFT was 
applied to the extended Hubbard model in seminal work 
by Sun et al. (37l . |49| . but only recently have numer- 
ical techniques for the solution of dynamical impurity 
models |44|-|46[ been sufficiently advanced to extract real- 
frequency information from such calculations [HI, l39j |. 
We employ the techniques of the latter two works (in 
particular a continuous-time quantum Monte Carlo im- 
purity solver within the hybridization expansion 
but implement them for the realistic Hamiltonian de- 
rived above. Moreover, we go beyond the 'standard' ex- 
tended Hubbard model and do not restrict the range of 
the non-local interaction terms. Rather, we include the 
entire 1/r tail by means of an Ewald-type lattice sum. 
In Fig. [5] we show momentum- resolved spectral functions 
from GW+DMFT for all compounds in our series: As ex- 
pected from the large onsite interactions compared to the 
bandwidth we obtain insulating spectra for all four com- 
pounds. Interestingly, however, for the Pb compound, in 
contrast to the other three systems, we find two stable 
solutions at the temperature of our study (T = 116K) - 
one metallic and one insulating. This indicates that we 
are in a coexistence region of a first order phase transition 
similar to that seen in the extended Hubbard modelpjfjj]. 

In all compounds the upper and lower Hubbard bands 
show substantial dispersion following the bare bandstruc- 
ture, as expected on general grounds. The insulating gap 
decreases within the series and we can estimate from the 
center of mass of the Hubbard bands values of: 1.3eV (C), 
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FIG. 3. (Color online) Left hand side: Frequency-dependent 
U(uj) (calculated from GW+DMFT) including both, insulat- 
ing and metallic cases for the Pb system. Right hand side: 
Imaginary part of the charge-charge susceptibility along the 
usual path in the Brillouin zone. 



0.8eV (Si), 0.7eV (Sn), and 0.5eV (Pb). However, specifi- 
cally for the Si(lll):{Sn,Pb} we find substantial spectral 
weight already at > — 0.2eV. Given this small gap, a siz- 
able temperature dependence can be expected. We have 
extracted the value of the local (i.e., k-integrated) spec- 
tral function at the Fermi-level [50] (see Fig.[2]bottom left 
panel). While for Si(lll):C the spectral weight transfer 
to the Fermi energy with temperature is negligible as ex- 
pected from the spectral function, specifically Si(lll):Si 
and most of all Si(lll):Sn display significant transfer of 
spectral weight on a temperature scale from 50K to room 
temperature 300K. 



Photoemission experiments for Si(lll):Sn jlfj, Il8 | 
(and, possibly [Hlj], for Si(lll):Pb 11|) observe, indeed, 



such a temperature dependence and agree well with our 
results, both, concerning the gap size and temperature 
scale. Our results - obtained without any adjustable pa- 
rameters - also stand as a theoretical prediction for more 
extensive studies on Si(lll):Pb and the (experimentally 
so far not studied) Si(lll):C compound. Next, we an- 
alyze the spectral functions in view of the interaction 
strengths calculated by cRPA (see Tab.HJ). The gap sizes 
no longer reflect the energy-scale of the onsite interac- 
tion Uq but are reduced due to non-local interactions 
which screen the local interaction by non-local charge 
fluctuations. This physics is naturally present in the 
GW+DMFT scheme, where non-local effects are incor- 
porated into an effective retarded onsite interaction IA{uj) 
(plotted in the left panel of Fig. [3]) . The shape of this 
quantity is reminiscent of screened interactions as calcu- 
lated, e.g., within the cRPA[32T]. where retardation ef- 
fects result from downfolding of high-energy degrees of 
freedom. The GW+DMFT U{lo) can be viewed as an 
effective interaction, where the dynamical character re- 
sults from downfolding non-local degrees of freedom into 
a local quantity. At large frequencies, screening is not ef- 
ficient and, hence, U{u) = oo) = Uq. On the other hand, 
the static value U(uo = 0) can be significantly reduced 
(up to nearly a factor of 2 for Si(lll):Pb). The latter 
sets the energy scale for the gaps we observe in the spec- 
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tral function. The transition between unscreened high- 
frequency behavior and the static value takes place at an 
energy scale ujq (plasmonic frequency) characteristic of 
the non-local charge fluctuations. The strikingly different 
behavior of the dynamical effective interactions M(uS) re- 
flects the observed materials trend: Si(lll):C [Si(lll):Si] 
is [nearly] unaffected by non-local interaction terms and 
there is barely any screening. The remaining two com- 
pounds show, however, large effects. The static values 
U (u> = 0) are reduced compared to the onsite interac- 
tion to 0.84cV for Si(lll):Sn and to 0.67cV (0.54cV) for 
the insulating (metallic) solution for Si(lll):Pb which 
leads to the reduced gap sizes. Moreover, plasmonic res- 
onances at energies between 0.6eV and 0.8eV stress the 
importance of non-local interactions/chargc-fluctuations 
for these systems. 

Besides leading to a retarded, frequency-dependent in- 
teraction, the non-local charge fluctuations signal ten- 
dencies towards a charge-ordered (CO) state. Analyzing 
the momentum-dependence of the imaginary part of the 
charge-charge response function Imx(k, u> = 0) for the 
high symmetry points of the Brillouin zone, shown in 
Fig. 02 we find for the different materials very particu- 
lar behavior. The local double occupancy, which corre- 
sponds to the integral of the plotted quantity over all 
momenta, becomes larger towards the end of the series. 
Most interesting is the case of metallic Si(lll):Pb for 
which we find a distinct structure within the Brillouin 
zone: The maximum of Imx(k, uj — 0) at the K sym- 
metry point indicates strong charge fluctuations of the 
so-called 3x3 symmetry, sketched in Fig, @] This order 
might eventually be frozen in to form a charge ordered 
ground state which is actually seen in scanning tunnel- 
ing microscopy for this material An insulating charge 
ordered ground state of 3 x 3 symmetry is, in fact also 
seen in Ge(lll):Sn(47| where a concomittant structural 
distortion (vertical displacement of adatoms) of the same 
symmetry is seen - our results show, that the instability 
in the correlated electronic response function is a good 
candidate for the key player of this feature. 

We can summarize our results by drawing a schematic 
phase-diagram as a function of the strength of local and 
non-local interactions (represented by the value of U\) 
as we show in Fig. |4] For zero non-local interactions 
our phase diagram describes the Mott-Hubbard metal to 
insulator transition. The adatom systems are placed at 
about 0.5eV of non-local interaction strength. However, 
due to the difference in the onsite term Uq their respective 
position in the phase diagram and, hence, their ground 
state character is different: Si(lll):C is deep in the Mott 
phase with a charge localization defined by one electron 
per adatom-site. 

The Si(lll):Si compound]!!! is also of Mott type with 
only small values for the double occupancy and little 
effect of plasmon excitations. However, Si(flf):Sn and 
most dramatic Si(lll):Pb (which is actually already in a 



> 

.Si. 

c 
o 

'J 

u 

CO 
>_ 

CD 

_E 

TO 
U 
O 
■ 

c 
o 
c 




Pb Si* C 

o o o o 

Sn 



V3x V3 

o o 
Mott o\p\o 
o o 



0.4 0.6 0.8 1.0 1.2 1.4 
local interaction [eV] 



FIG. 4. (Color online) Schematic local/nonlocal- interaction 
phase-diagram: The black-bordered circles mark the positions 
of the adatom systems of our study. Straight lines are guide 
to the eyes and blurry color indicates coexistence regions. 
Within the localized (i.e. insulating) phases CO and Mott, 
small sketches indicate the shape of the surface unit-cell. 



coexistence region) are much closer to a phase boundary 
to a metallic phase. Even more peculiar is the obvious 
tendency of Si(lll):Pb towards a charge- ordered phase 
of 3 x 3 symmetry indicated by the white region in our 
phase diagram. 

In conclusion, we have set up a fully self-consistent 
GW+DMFT scheme for the realistic treatment of corre- 
lated surface systems to address the electronic properties 
of the a-phases of adatoms on the Si(lll) surface. We 
reported on the ab initio construction of the materials- 
specific low-energy Hamiltonians and, most importantly, 
on the respective interaction parameters including the 
long-range Coulomb tail. From these it becomes clear 
that for the adatom systems taking into account non- 
local interaction effects is mandatory. We have solved 
the derived many body Hamiltonians and discussed our 
finding for momentum-resolved spectral functions, to be 
compared to ARPES spectra. Without adjustable pa- 
rameters we reproduced available experimental findings, 
or, in (most) cases where experiments are missing, made 
predictions. Specifically, the ARPES spectra for the se- 
ries, as well as the charge order instabilities in the case 
of Si(III):Pb are key conclusions/predictions which can 
provide guidance for further experimental and theoretical 
studies of semiconductor adatom structures. 
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